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Quantum optimal control theory is applied to two and three coupled Josephson charge qubits. It is 
shown that by using shaped pulses a cnot gate can be obtained with a trace fidelity > 0.99999 for the 
two qubits, and even when including higher charge states, the leakage is below 1%. Yet, the required 
time is only a fifth of the pioneering experiment 1] for otherwise identical parameters. The controls 
have palindromic smooth time courses represent able by superpositions of a few harmonics. We 
outline schemes to generate these shaped pulses such as simple network synthesis. The approach is 
easy to generalise to larger systems as shown by a fast realisation of Toffoli's gate in three linearly 
coupled charge qubits. Thus it is to be anticipated that this method will find wide application 
in coherent quantum control of systems with finite degrees of freedom whose dynamics are Lie- 
algebraically closed. 

PACS numbers: 85.25.Cp, 82.65.Jn, 03.67.Lx, 85.35.Gv 



In view of Hamiltonian simulation and quantum com- 
putation recent years have seen an increasing amount 
of quantum systems that can be coherently controlled. 
Next to natural microscopic quantum systems, a par- 
ticular attractive candidate for scalable setups are su- 
perconducting devices based on Josephson junctions 0- 
Due to the ubiquitous bath degrees of freedom in the 
solid-state environment, the time over which quantum 
coherence can be maintained remains limited, although 
significant progress has been achieved 0,3- Yet, it is a 
challenge how to produce accurate quantum gates, and 
how to minimize their duration such that the number 
of possible operations within T2 meets the error correc- 
tion threshold. Concomitantly, progress has been made 
in applying optimal control techniques to steer quantum 
systems 5j in a robust, relaxation- minimising 6] or time- 
optimal way 7] . Spin systems are a particularly powerful 
paradigm of quantum systems 8] : under mild conditions 
they are fully controllable, i.e., local and universal quan- 
tum gates can be implemented. In N spins- ^ it suffices 
that (z) all spins can be addressed selectively by r/-pulses 
and {it) that the spins form an arbitrary connected graph 
of weak coupling interactions. The optimal control tech- 
niques of spin systems can be extended to pseudo-spin 
systems, such as charge or flux states in superconducting 
setups, provided their Hamiltonian dynamics can be ap- 
proximated to sufficient accuracy by a closed Lie algebra, 
e.^., in a system of N qubits 0u(2^). 

As a practically relevant and illustrative example, we 
consider two capacitively coupled charge qubits con- 
trolled by DC pulses as in Ref. The infinite- 



dimensional Hilbert space of charge states in the device 
can be projected to its low-energy part defined by zero or 
one excess charge on the respective islands 2]. Identify- 
ing these charges as pseudo-spin states, the Hamiltonian 
can be written as i^tot = ^drift + ^controb where the 
drift or static part reads (for the constants see caption to 
Fig. 1) 
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while the controls can be cast into 
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-ff control = (—"•92+ E^iUgi ] (ay ® 1) 
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The control amplitudes Ug 



1,2 are gate charges 



controlled by external voltages via Ugi, = VguCgu/2e. 
They are taken to be piece-wise constant in each time 
interval tk. This pseudo-spin Hamiltonian motivated by 
Ref. also applies to other systems such as double quan- 
tum dots ^ and Josephson flux qubits 10], although in 
the latter case the controls are typically r/-pulses. 

In a time interval tk the system thus evolves under 

(k) (k) 

H^^( = i^drift + ^control- ^^^^ ^ scqucuce 

of control amplitudes for the times ti, ^2, . . . , t/^, . . . , ^at 
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FIG. 1: (Color online) Fastest charge level controls obtained 
for realising a CNOT-gate on a pair of coupled charge qubits 
(left part: control qubit, right part: working qubit). The 
total gate charges for the qubits are Ug^ — n^^J + Sugu with 
= 1, 2. Here, n^g-^ — 0.24, rSg2 =0.26 and the qubit energies 
Eci/h = 140.2 GHz, Ec2/h = 162.2 GHz, Eji/h = 10.9 GHz, 
Ej2/h — 9.9 GHz, and Em/h — 23.0 GHz were chosen in 
accordance with the experimental values given in 1]. The 
50 piecewise constant controls are shown as bars; the trace 
fidelity is ^ |tr{[//^rgetf^T}| > 1 - 10"^ Red lines give the 
analytic curves in Eqn. El while the blue ones superimposed 
show a pulse actually synthesised by an LCR-filter (see a later 
section and Fig. 3). 

such as to maximise a quality function, here the overlap 
with the desired quantum gate or element of an algo- 
rithm /7target- Morcovcr, for the decomposition of Ut = 
^-UnHn^-Un-iHn-i . . . ^-itkHk . . . ^-itiH^ ^^^^ available 

controls {H^} to be timeoptimal, T := Yl'^=i has to be 
minimal. The gate fidelity is unity, if ||/7t — tA;arget|l2 = 
= \\Ut\\1 + ||t/target|l2 " 2RetT{uL^M. Maximis- 
ing Re tr{/7/g^rgg^/7T} for fixed T can readily be solved by 
optimal control: Let h{U{tk)) := Re tr{A'''(t/e)(— i(i^d + 
^UjjHjj))U{tk)} with the Lagrange-type adjoint system 
\(t) following the equation of motion \{t) = —i{H(^ + 
^UjjHjj)\{t). Pontryagin's maximum principle requires 
dh/duiy = Re tT{X'^ {—iHj^)U} = thus allowing to im- 
plement a gradient-flow based recursion. For the am- 
plitude of the control in iteration r + 1 at time 
interval tk one finds with e as a suitably chosen step 
size n^Qj^^\tk) = rihv itk) + (^fc) derived in 

Refs. [lllll^- Here T is the shortest time allowing for a 
given fidelity numerically. 

We now turn to the discussion of our numerical re- 
sults. We have used parameter values from the exper- 
iment 0. Variation of these values should change de- 
tails of the result, but not its overall structure. Fig. 1 
shows the fastest decompositions obtained by numeri- 
cal optimal control for the CNOT gate into evolutions 
under available controls (Eqns. 1 and 2). In contrast 
to the 250 ps in Ref. 0, T = 55 ps suffice to get 
\\Ut — t^target||2 = 5.3464 X 10~^ Corresponding to a tracc 

fidelity of jrHuL^Ml > 1 - 10-'- 

Beyond the efficient and accurate implementation, this 



result provides physical insight: our pulse essentially ac- 
comodates all terms of the standard CNOT pulse sequence 
for this coupling 2] such that different terms in the total 
Hamiltonian act in parallel instead of sequentially. For a 
CNOT, the duration T = 55 ps has to accomodate at least 
a ^ rotation under the coupling Hamiltonian {^(Jz^ctz) 
lasting 21.7 ps concomitant to two ^ x-rotations under 

the second of the drift components (|cri^^ with u = 1,2) 
requiring 22.9 and 25.3 ps, respectively. Thus, unlike in 
NMR, the time scales of local and non-local interactions 
are comparable. Assume in a limiting simplification that 
two ^ x-pulses are required, the total length cannot be 
shorter than 50.6 ps. Our solution is close to this infi- 
mum. Note that a duration of T = 55 ps also implies that 
the trajectory of the coherent evolution does not have to 
be a geodesic in the Weyl chamber (compare ref. 7]), as 
shown in the supplement Fig. 5. Moreover, the evolution 
times for single components do not add up, indicating 
parallel evolution of different interactions. 

The supplementary material illustrates how the se- 
quence of controls (Fig. 1) acts in a quasi-continuous way 
on specific input states: Suppl. Fig. 1 gives the evolu- 
tion of a product state, |0) = |0)|0). The representation 
of the reduced states in their local Bloch spheres shows 
how the control qubit undergoes a closed loop, while 
the working qubit is inverted as expected. As demon- 
strated in Suppl. Fig. 2, a maximally entangled Bell state 
1?/^+) = -^(|0)|0) + |1)|1) evolves from the centre of the 
Bloch spheres (indicating maximal entanglement) into a 
product state on the surface of the sphere. 

Note that the time course of controls in charge 
qubits turns out palindromic (Fig. 1). Self-inverse 
gates {Ugg^^Q = 1) relate to the more general time-and- 
phase-reversal symmetry (TPR) observed in the con- 
trol of spin systems [l^: for example, any sequence 
g-*t^cr^g-ztyayg-2t^cr^ -g inverted by transposition con- 
comitant to time reversal tiy —tiy and ay ~<^y 
Since the Hamiltonians in Eqns. 1-2 are real and sym- 
metric, they will give the same propagator, no matter 
whether read forward or backward. 

The pulse is not very complicated. Interestingly, the 
time-course of the controls on either qubit (i/ = 1, 2) can 
be written as a sum of 6(7) harmonic functions 

5(6) ^ 
^^^W = X^^^O') cos (27rcj^(j)- + 0^(j)) . (3) 

j=0 

The constants from Tab. 1 in the supplementary material 
give a high accuracy (x^ = 0.008231; 0.003668 for the 
channels 1 and 2, respectively). 

This representation reflects the simplicity and the 
modest bandwidth of the pulses obtained. The low band- 
width allows to maintain a high fidelity even if leakage 
levels formed from higher charge states of the qubit sys- 
tem are taken into account: we now explicitly apply the 
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FIG. 2: (Color online) Spectroscopic explanation of the high quality of the control sequences of Fig. 1: the spectral overlap 
of the Fourier- transforms (right walls) of the controls of Fig. 1 with the energy differences corresponding to the one-charge 
transitions into leakage levels (solid lines on the surface) is small at charges around Ugu — 0.2 with ly — 1,2. Intensities 
at allowed (solid lines) vs forbidden transitions (broken lines) into leakage levels are given in terms of the transition-matrix 
elements f\Hc^i) normalized by the charging energies Ed {Ec2) in the 3D representation: the working transitions (blue) are 
of the same probability as the allowed transitions into leakage levels (red) that have no overlap with the excitation bandwidth 
of the pulses, while the forbidden ones are too weak to show up at all. (details in the text). 



full pulse to the extended Hamiltonian obtained by map- 
ping the full Hamiltonian [ij to the subspaces of — 1, . . . , 2 
extra charges per island. The two-qubit cnot gate is 
then embedded into the group SU{1&). Even then the 
propagator generated by the above controls projects well 
onto the CNOT gate giving a trace fidelity > 0.99. The 
good result may be astounding at first sight, however, 
it can be understood by relating the limited bandwidth 
to the matrix elements, which both control the transi- 
tion rate: while the one- charge transitions to the leakage 
levels like | — 1) ^ |0) and |2) ^ |1) are allowed by 
the Josephson coupling, the two- charge transitions like 
I — 1) ^ |1) and 1 2) ^ |0) are forbidden in terms of the 
transition probabilities f\Hc^i) as can be seen from 
Fig. 2. Moreover, note that the charge levels of Fig. 1 are 
mostly around = 0.2 thus contributing to the working 
transition |0) ^ while the 'spectral overlap' of the 
Fourier-transform of the time course of the controls with 
energy differences corresponding to the potentially dele- 
terious one-charge transitions in Fig. 2 is small. Hence 
there are simple spectroscopic arguments for the high fi- 
delity obtained by the controls of Fig. 1. 

Furthermore, even the time courses starting out with 
any of the four canonical two-qubit basis vectors hardly 
ever leave the state space of the working qubits: at any 
time the projections onto the leakage space do not exceed 
0.6 %. Choosing initial states from the Bell basis entails 
even less leakage. Note that explicitly taking into account 
the leakage levels during optimisation is expected to im- 
prove the quality even further. Thus, also pseudospin 
systems like ours, which involve a low-energy projection 
disregarding leakage levels, can be controlled with high 
accuracy and quantum computing is not strictly limited 



to coupled two- level systems. 

Generating these pulse shapes experimentally is a chal- 
lenging but possible task. Note that the length of the 
pulse is given by the coupling strength as discussed above 
and hence can be extended by lowering the coupling. 

In the pertinent time scale, there are no devices comer- 
cially available for generating arbitrary wave forms with 
the same capabilities as NMR-spectrometers. High-end 
commercial pulse generators as well as custom-built ones 
are close to the necessary specifications 0, . Pulses 
can be formed by superimposing short pulses of shapes 
easy to generate with different heights, widths, and de- 
lays. The two main candidates for this approach are (z) 
Gaussian pulses |2| , which can be generated at room tem- 
perature and which run nearly undistorted through the 
necessary cryogenic filtering and (n) SFQ pulses, which 
can be generated on chip (hence avoiding the filters) us- 
ing ultrafast classical Josephson electronics 0, ^3 • 

We would like to exemplify a well established tech- 
nique, shaping in Laplace space, to generate these pulses. 
The idea resembles the approach of femtosecond quan- 
tum chemistry: we start with an input current pulse 
/in(t) shorter than the desired one of a shape which is 
arbitrary as long as it contains enough spectral weight 
at the harmonics necessary for the desired pulse. Such 
pulses are readily generated optically or electrically and 
have, without shaping, already been applied under cryo- 
genic conditions 15]. This pulse is sent through a discrete 
electrical four-pole, whose transfer function Z12 is de- 
signed such that the desired pulse is found at the output. 
We have carried out this idea for a rectangular pulse of 
length Tr = l.lps as an input and our two gate pulses as 
output. We have developed a transfer function in Laplace 
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FIG. 3: (Color online) Characterization of a filter shaping 
the pulse on the second gate. The bars show the poles Si of 
the transfer function in the Laplace plane. Poles outside the 
negative imaginary axis also lead to the complex conjugate 
pole and can be implemented by an LCR-Filter. The height 
of the bars show the modulus of the residue in this pole. 



space Zi2{s) by fitting Vg{s) = Zi2{s)Iin{s) (see Fig. 3). 
Owing to causality, the poles of Z12 are either on the 
negative real axis or in conjugate pairs of poles on the 
left half plane. Each conjugate pair corresponds to an 
LCR-filter stage whereas each real pole corresponds to 
an RC lowpass- filter. It turns out, that good agreement 
can be achieved with 8 LCR filters and two low-pass fil- 
ters, following the standard rules of circuit synthesis 

The pulses are very close to the desired ones, see Fig. 1, 
and a trace fidelity of 94 % can be achieved for the entire 
CNOT. Clearly, the quality can be further improved with 
more refined technology. 

The filter as well as the pulse design are ready to acco- 
modate the experimental necessities. On the one hand, 
due to unavoidable fabrication uncertainties, the opti- 
mum pulse will look slightly different for each individual 
pair of qubits. Realistically, the matrix elements of the 
total Hamiltonian Eqs. Q first have to be deter- 

mined spectroscopically, then our algorithm has to be 
run to find the optimum pulse shape. This is done on 
a regular PC in ca. 30 seconds. Secondly, for adjusting 
the filtering circuit which can be put at room tempera- 
ture, one has to take into account the transfer function 
through the filters of the cryostat and to the sample. 
This contribution ^sample to the total transfer function 
will most realistically be measured using a capacitor of 
the same geometrical dimensions of the qubits as a probe. 
As long as this does not block the relevant frequencies, 
i.e., if the setup has sufficient bandwidth, ^sample can be 
accounted for when adjusting an additional pulse shap- 
ing filter such that the total transfer function shapes the 
correct pulse. Note, that our method also applies to con- 
trol by microwave Rabi-type pulses, where pulse shaping 
appears to be easier as time scales are usually longer. 
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FIG. 4: Fastest charge level controls obtained for realis- 
ing a Toffoli-gate on a linear chain of charge qubits cou- 
pled by nearest-neighbour interactions. The piecewise con- 
stant controls are shown as bars. The trace fidelity is 
|tr{[//g^rgg^[/T}| > 1 — 10~^. Here, the parameters are 
Eci/h = 140.2 GHz, Ec2/h = 120.9 GHz, Ecs/h = 184.3 
GHz, Eji/h = 10.9 GHz, Ej2 = /h9.9 GHz, Ejs/h = 9.4 
GRz,Emi,m2/h = 23 GHz, n°,i = 0.24, n°,2 = 0.26, n^,3 = 
0.28. 



Likewise, in a system of three linearly coupled charge 
qubits, one may decompose the TOFFOLI gate into ex- 
perimentally available controls. 

This result highlights that due to the comparatively 
strong qubit-qubit interactions in multiqubit setups, the 
direct generation of three-qubit gates is much quicker 
than its compostion into elementary universal gates , e.g. 
decomposing a Toffoli into 9 CNOTS in a linear spin 
chain: the speed-up is by a factor of 2.8 compared to 9 of 
our CNOTS and by a factor of 13 compared to Nakamura's 
CNOTS 1]. 

This also holds when developing simple algorithms 
on superconducting qubit setups: a minimization algo- 
rithm for searching control amplitudes in coupled Cooper 
pair boxes has been applied in 20], however, in that ap- 
proach, the numerical optimization was restricted to only 
a few values. In Ref. |2l|, a pulse sequence generating a 
CNOT with fixed couplings has been invented, which uses 
hard RF pulses instead of our shaped pulse and turns out 
to be much longer, thus leads to serious conflicts with de- 
coherence. 

In conclusion, we have constructed pulses for the real- 
ization of fast high-fidelity quantum logic gates in super- 
conducting charge qubits. The optimum pulses are al- 
ways palindromic, owing to the time-reversal invar iance 
of these pseudo-spin Hamiltonians. The simplicity of the 
pulse shape results in low bandwidth and thus low leak- 
age to higher states and the setup necessary to generate 
such pulses is of modest complexity. 

We are indebted to Navin Khaneja for continu- 
ous stimulating scientific exchange. We thank M. 
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DETAILED PRESENTATION OF THE PULSE AND THE QUBIT DYNAMICS 



In this supplement, we provide more insights into the numericahy obtained pulse and the resulting qubit dynamics 
as an additional visualization. Table I contains the precise numerical settings for the Fourier decomposition of the 
pulse shown in Fig. 1 and Eq. 3 of the paper. Figure Q shows the qubit dynamics on the reduced Bloch-spheres of 
the two qubits, i.e. the left Bloch spere, belonging to the control qubit, qubit 1, shows the spin projections of the 
reduced density matrix, pi = Tr2/?, and vice versa. Even if the two-qubit state is pure, the reduced Bloch vector 
may be inside the sphere, which hints at entanglement between the qubits, as seen in Fig. [21 The experimental pulse 
shown in Fig. 01 requires a much longer trajectory on the sphere. The representation in the Weyl chamber. Fig. [^1 
provide a complementary visualization. It represents the generic, irreducible two qubit part alone, i.e. each point in 
the chamber is invariant under single-qubit rotations. Details of this representation can be found in Ref. 0. 

PULSE SHAPING HARDWARE 

This section details the pulse shaping scheme outlined in the main manuscript. The data in Fig. 3 of the paper 
and in Fig. of this supplement have been obtained as follows: We have fitted a rational function ^12(5) in Laplace 
space, such that Vout{s) = Zi2{s)Iin{s) where /in is a Ips current pulse and Ug^i = CG,i'K)ut,i/2e for the two qubits. 






FIG. 1: Evolution of the product state |O(0)) = |0)|0) under the optimised controls resulting in |G(T)) = |0)|1). The evolution 
<t <T with T = 55 ps is represented by the reduced states trB\0{t)){0{t)\ (left sphere) and trA\0{t)){0{t)\ (right sphere) 
on the respective local Bloch spheres. The blow-up on the left shows the top of the left Bloch sphere. 




FIG. 2: Evolution of the Bell state |$+) = ^ (|00) + |11)) into the final state ^ (|01) + |11)) (filled red dots). The Bell state 
is maximally entangled and hence has local representations in the centre of the respective Bloch spheres, while the final state 
is a product state represented by points (filled red dots) on their surfaces. The projection on the left is a view from the top 
onto the plane inserted into the left Bloch sphere. 
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TABLE I: Parameters giving the envelope to the control amplitudes Uguit) for the two qubits = 1, 2 as in Eqn. 3 in the main 
text. T = 55 ps. The goodness of fit is determined by = 0.008231 for qubit 1 and = 0.003668 for qubit 2. 
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i = 1, 2. This function is represented best by its residue decompostion Z12 = Xli ^^7"- With this decomposition, there 
are a number of approaches to design a lumped circuit with this transfer function, such as the method of Gewertz 
1^. This methods systematically eliminates poles and introduces loops in the electrical circuit: An LCR-loop for each 
pair of complex conjugate poles, an RC-filter for each pole on the real axis. Thus, the degree of the polynomial in 
the denominator gives a clear view on the size of the necessary circuit. In reality, the transfer function from the pulse 
shaping circuit, which can conveniently be placed at room temperature, to the sample is not smooth. By using a passive 
classical pickup element simulating the qubit, e.g. a capacitor in the charge qubit case |1||. This transfer function 




FIG. 3: For comparison: same evolution of the product state |O(0)) = |0)|0) as in Fig. 1, but using the pulse of the experiment 
Ref. Q. The evolution < t < T with T = 255 ps is represented by the reduced states trs|e(t))(e(t)| (left) and trA|B(t))(e(t)| 
(right) on the respective local Bloch spheres. The trajectory completes two full circles (see inset) before reaching its final state 
near the south pole. 






FIG. 4: Evolution of the Bell state |$+) = ^ (|00) + |11)) as Fig. 2, but using the control of the experiment 1]. Parameters 
of that pulse on the second qubit are: 5ng2 — 0.25, 255 ps total length, with 40 ps rise time and 40 ps fall time, digitisation: 
1000 points. Note the different final states as compared with Figs. 1 and 2 indicative of a resulting gate whose matrix elements 
coincide with the proper cnot in absolute value, but not in phase. (Actually, the phase deviations are not uniform throughout 
the elements). 
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(Jt,0,0) 

FIG. 5: Trajectory for the optimised cnot in the Weyl chamber (compare ref. 0). Starting at the origin, the trajectory is a 
non-geodesic smooth curve (see explanation in the main text) which ends at the point (-1, 0, 0) as expected for cnot. 




FIG. 6: Pole structure of the transfer functions necessary for shaping the pulses on both the control (left) and target (right) 
qubit. Parameters correspond to Fig. 3 of the main paper. 

can be measured. In the linear case, it can be expressed as another four-pole impedance matrix Zgampie- The total 
transfer function of the series configuration of those four-poles will be Zyi = ^12, sampie^i2,fiiter/(^22, filter + ^ii, sample)- 
This outlines the statement in the text, that unless the transfer function to the sample is not filtering out the relevant 
frequencies (i.e. becomes small for values of s important to V^ut), it will be possible to design an appropriate filter 
taking into account the properties of the experimental setup. 
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